Magnetization relaxation in (Ga,Mn)As ferromagnetic semiconductors 
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We describe a theory of Mn local-moment magnetization relaxation due to p-d kinetic-exchange 
coupling with the itinerant-spin subsystem in the ferromagnetic semiconductor (Ga,Mn)As alloy. 
The theoretical Gilbert damping coefficient implied by this mechanism is calculated as a function of 
Mn moment density, hole concentration, and quasiparticle lifetime. Comparison with experimental 
ferromagnetic resonance data suggests that in annealed strongly metallic samples, p-d coupling con- 
tributes significantly to the damping rate of the magnetization precession at low temperatures. By 
combining the theoretical Gilbert coefficient with the values of the magnetic anisotropy energy, we es- 
timate that the typical critical current for spin-transfer magnetization switching in all-semiconductor 
trilayer devices can be as low as ~ f0 5 A cm -2 . 

PACS numbers: 73.20.Mf, 73.40.-c,85.75.-d 



I. INTRODUCTION 

The Gilbert coefficient describes the damping of small- 
angle magnetization precession and is one of the key pa- 
rameters that characterizes collective magnetization dy- 
namics in a ferromagnet. Early theories of magnetization 
dynamics in transition metals viewed exchange coupling 
(oc S • s) between local moment d-shell spins S and itin- 
erant s-p band spins s as a key relaxation mechanism. 1 
We now recognize that this model needs to be elaborated 
for transition metals to account for the itinerant charac- 
ter of their d-electrons. Models of d-shell local moments 
that are exchange-coupled to itinerant s-p band electrons 
have, however, been resurrected recently because they 
provide a good description of ferromagnetism in many 
diluted magnetic semiconductors (DMSs), (Ga,Mn)As 
in particular^' 3 Exchange-coupling between local mo- 
ments and itinerant electrons should also contribute sig- 
nificantly to Gilbert damping in these new ferromagnetic 
systems. The elementary process for this damping mech- 
anism is one in which a local-moment magnon is annihi- 
lated by exchange interaction with a band electron that 
suffers a spin-flip. This process cannot by itself change 
the total magnetic moment since the exchange Hamilto- 
nian commutes with the total spin S + s. Net relaxation 
of the magnetization requires another independent pro- 
cess in which the itinerant electron spin relaxes through 
spin-orbit interactions. 

Recent experiments and ab-initio calculations* have es- 
tablished that for doping levels up to several per cent, 
a substitutional Mn impurity in GaAs introduces five 
strongly localized d-electrons and a delocalized hole in 
the As/Ga p-band, and that ferromagnetic coupling be- 
tween the S — 5/2 Mn moments is mediated by the 
p-d kinetic-exchange. Hence, we model* 5 - the electronic 
structure of the free carriers by that of the host mate- 
rial, implicitly assuming a shallow acceptor picture. The 



free carrier quasiparticles are p-d exchange coupled to 
the local moments with strength 3 J p d ~ 55 meV nm 3 , 
and have a finite life time that can be estimated pertur- 
batively. This theoretical picture leads to an accurate 
description of many thermodynamic and transport prop- 
erties of optimally annealed (Ga,Mn)As samples, such 
as the measured transition temperature?^ the anoma- 
lous Hall effect ?M anisotropic magnetoresistancejS and 
magneto-optical properties'**^ Particularly important in 
justifying the present theory are results for the magneto- 
crystalline anisotropy?A2, spin-stiffnessjii* 1 ^ and Bloch do- 
mam width 13 that all agree well with experiment. These 
parameters follow from the long-wavelength limit of the 
theory of the Mn spin-wave dispersion, and reflect the 
retarded and non-local character of the valence-band- 
hole mediated interactions between Mn moments^* The 
Gilbert damping of magnetization precession discussed 
here is the aspect of this long-wavelength collective mag- 
netization dynamics that is most directly dependent on 
valence-band spin-orbit coupling. 

In Section II of this paper, we present a fully mi- 
croscopic theory of the kinetic- exchange contribution to 
the Gilbert coefficient in DMSs. By comparing the lin- 
ear response predicted by the classical phenomenologi- 
cal Landau-Lifshitz-Gilbcrt (LLG) equation with micro- 
scopic linear response theory, we identify the Gilbert co- 
efficient with the dissipative part of a susceptibility dia- 
gram. 

In Section III we present experimental ferromagnetic 
resonance (FMR) data— ill recorded as a function of tem- 
perature and external magnetic field strength and ori- 
entation. Since the frequency dependence of the FMR 
linewidth is not available fi 2 * 1 * we are unable to ex- 
perimentally decouple inhomogeneous FMR broadening 
and the intrinsic Gilbert damping contributions to the 
linewidth to make a quantitative comparison with the 
theory. Nevertheless, the data indicate that the magnetic 
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inhomogeneity contribution is largely suppressed in the 
more metallic, annealed samples and that much of the 
observed low-temperature FMR linewidth in these sam- 
ples can be explained by damping of the magnetization 
precession mediated through the p-d coupling. 

By adding a spin-torque termi£ to the LLG equation, 
we estimate in Section IV that the typical critical current 
for magnetization switching due to spin-transfer torques 
in an all-semiconductor trilayer device consisting of mag- 
netically "soft" and "hard" DMS layers separated by a 
non-magnetic spacer will be ~ 10 5 Acm~ 2 . In metals, 
this spin-transfer effect is currently the focus of a consid- 
erable experimental^ and theoretical research. Spin- 
transfer switching has not yet been demonstrated in all- 
semiconductor systems, but the effect promises to have a 
richer phenomenology in this case because of the flexibil- 
ity of semiconductor ferromagnet materials, and because 
of the possibility of combining spin-transfer with other 
semiconductor circuit functionalities. 1819 The relatively 
low critical currents we predict for semiconductors may 
also circumvent the incomplete magnetization switching 
encountered in metallic magnetic tunnel junctions that 
occurs due to the interference of strong self-field effects 
with the spin-transfer torques^ 

The paper concludes with a brief summary of our re- 
sults. 



II. THEORY OF THE GILBERT DAMPING 

Semiclassical LLG Linear Response. The phenomeno- 
logical LLG equation for collective magnetization dynam- 
ics is 
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where M is the local Mn moment magnetization, B e // = 
—dE/dM is the effective magnetic field, g is the Lande 
g- factor, /is is the Bohr magneton, and a is the phe- 
nomenological Gilbert damping coefficient. Unless a de- 
pends strongly on the orientation of the magnetization 21 
or if the magnetization is not fully aligned with the ex- 
ternal static magnetic field, the Gilbert damping rate, 
observed in experiment through a frequency-dependent 
FMR linewidth, is independent of the static field and 
of the details of magnetic anisotropics present in the 
sample. 22 This allows us to assume in this section a sim- 
ple geometry in which the anisotropy fields are repre- 
sented by a single, uniaxial anisotropy energy density 
coefficient U. For small fluctuations of the Mn magneti- 
zation orientation around the easy axis, Eq. Q can be 
used to derive a phenomenological response function of 
the magnetic system to a weak transverse field. For zero 
external static magnetic field the corresponding inverse 
susceptibility reads: 
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where U — U/(hNM n S), Nmv, — ^ x / a f c * s the density of 
uniformly distributed Mn moments in Gai-^Mn^As (a/ c 
is the GaAs lattice constant), and lu is the frequency of 
the external rf field perturbation. 

Microscopic theory. We derive the zero-temperature 
quantum response function from our effective Hamilto- 
nian theory and obtain a microscopic expression for a by 
equating the quantum mechanical response function to 
the classical one in the uniform u> — > limit. We start 
by writing a quantum analog of Eq. Q using the linear 
response theory, 

{M x {v,t))=- % -J dt' J dr , (^([M x (r,t),-M x (r',t') 
B x {t')]) + ([M x (r, t),-M v (r', t')B v (t')})) 6{t - t') , (3) 



which leads to the retarded transverse susceptibility: 

(4) 

Here i = x,y and 5j(r, t) = Mi(r,t)/(gfig) are the Mn 
transverse spin-density operators. 

To evaluate the correlation function J3J we use the 
Holstein-Primakoff boson representation of the spin oper- 
ators assuming small fluctuations around the mean-field 
ordered state, S+ = by/2N Mn S and S~ = b^2N Mn S 
(S x = (S+ + S-)/2, S y = (S+ ~ S-)/2i), and choosing 
the z-direction as the quantization axis. After integrat- 
ing out the itinerant carrier degrees of freedom within 
the coherent-state path-integral formalism of the many- 
body problem we obtain the partition function, Z = 
J T>[zz]exp(—S[zz]), with the action given to quadratic 
order in z and z (the complex numbers representing the 
bosonic degrees of freedom) by 



S[zz] = l//3V^2z(k,n m )(-iQ r , 

rn , k 



(5) 

In Eq. (JSJ), the first term is the standard Berry's phase 
and the second term is the itinerant carrier spin polariza- 
tion diagram for the high symmetry case where the cross 
terms of the form zz and zz vanish in S , [S^]iai 
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Here Q a (z,li) is the single-particle band Green's func- 
tion and <^ (k) are the band eigenstates. From the par- 
tition function we compute directly the imaginary time 
response functions at finite Matsubara frequencies which 
after their corresponding analytic continuations yield 
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Here II rGt (calculated below) describes mathematically 
the retarded interaction between the Mn bosonic degrees 
of freedom due to the p-d kinetic exchange with valence 
band holes. 

Connecting classical phenomenology and microscopic 
theory. Inverting the retarded susceptibility in Eq. (JJJ 
for the uniform (k = 0) precession mode, we obtain 
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Comparing Eqs. |JS} and (J2J, we obtain the micro- 
scopic contribution to the Gilbert coefhcient from kinetic- 
exchange coupling : 
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Note that x R y is explicitly equal to the off-diagonal 
element of x" 1 hi Eq. [2] and that, also consistently, the 
real part of n rct (cj), in the limit of co — > 0, gives the 
magnetocrystalline contribution to the anisotropy energy 
U, as explained in detail in Ref. 0. 

To compute II ret (a;), we take into account the effects 
of disorder present in the system perturbatively by ac- 
counting for the finite-lifetime of band quasiparticlcs, 
which for simplicity we characterize by a single num- 
ber r. The quasiparticle broadening T was chosen to 
be in the range estimated in previous detailed studies of 
transport properties of these systems, which achieve good 
agreement with experiment* 8 " The single particle Green's 
function for the valence-band quasiparticles is thus writ- 
ten as G a (k,z) = f^° oo dtj'/(2Tr)A a (ij',k)/(z - w') with 
a spectral function A a (e,k) = T/[(e - e aik ) 2 + T 2 /4]. 
In the present case we take the valence-band electronic 
structure to be described by the six-band Kohn-Luttinger 
Hamiltonian in the presence of an effective kinetic- 
exchange field h e ff = J p dNMn{S)*^ In a collinear fer- 
romagnetic state and zero temperature (S) = S, and we 
obtain 
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In choosing a spin- and band-independent lifetime, 
we are implicitly appealing to the dominance of spin- 
independent Coulomb scattering off Mn acceptors and in- 
tcrstitials as the dominant 8 scattering mechanism. Spin- 
orbit interactions enter through their presence in the in- 
trinsic bands rather than through spin-flip quasiparticle 



scattering events. In this model, the Gilbert coefficient 
in Eq. 1)1 0|l has intra-band and inter-band contributions 
that have qualitatively different disorder dependences, 
as illustrated in Fig. ^ Note that the terminology we 
use here recognizes that no band has spin-character suf- 
ficiently definite to justify the usual distinction between 
majority and minority spin bands. The intra-band term 
we refer to here would be spin-flip scattering within a 
given spin-split band in the more usual language and 
is present only because of intrinsic spin-orbit coupling 
in the host semiconductor bands. The intra-band con- 
tribution to a is proportional to 1/r at small T, i.e., 
proportional to the conductivity rather than the resis- 
tivity, and would dominate the damping if disorder were 
weak. The inter-band contribution to a, on the other 
hand, requires disorder to breach the wavevector separa- 
tion between different bands at the Fermi energy and is 
an increasing function of T. The overall dependence of 
the Gilbert coefficient on the sample's mobility is non- 
monotonic, as illustrated in Fig^ with the position of 
the minimum depending on both hole and Mn-moment 
densities, and on other parameters of the DMS system. 
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FIG. 1: Total Gilbert damping coefficient a, interband contri- 
bution «mtcr, and intraband contribution «i n tra as a function 
quasiparticle life-time broadening T for a carrier density of 
0.5nm" 3 and x = 8%. 

The complexity and tunability of the Gilbert coeffi- 
cient in DMSs is illustrated in Fig. |2 where we plot a as 
a function of the hole density for T — 150 and 50 meV 
and for Mn doping x =2 - 8%. These parameter values 
bracket the range typical for metallic (Ga,Mn)As DMSs. 
The theory curves in Fig. [21 predict that a increases with 
increasing hole density because of the higher quasiparti- 
cle density of states at larger densities. The dependence 
of the Gilbert coefficient on x is more complex. The 
prefactor h e ff in Eq. i|10|) reflects the proportionality of 
the kinetic- exchange coupling to the Mn spin density and 
causes the Gilbert damping implied by this mechanism 
to decrease with decreasing x at low Mn doping. This 



4 




Hole density (nm 3 ) Hole density (nm 3 ) 



behavior is clearly seen in Fig.^a) for x < 6 %. On the 
other hand, an opposite trend is predicted for higher Mn- 
moment densities where the effect of h e ff on the intra- 
and inter-band matrix elements in Eq. i|10[l takes over. In 
that case, larger h e ff values lead to a reduced spin mix- 
ing in the quasiparticle bands and, therefore, to smaller 
Gilbert damping rates. This implicit dependence of a on 
heff is more dramatic in higher-quality samples. 

We expect that the above kinetic- exchange mechanism 
of the Gilbert damping will dominate at low tempera- 
tures where other mechanisms such as magnon-magnon 
interactions vanish. At temperatures close to the Curie 
temperature, on the other hand, the contribution to 
magnetization precession damping due to the kinetic- 
exchange coupling can be ignored. The argument is based 
on an approximation that combines our zero-temperature 
microscopic theory of a with a finite-temperature mean- 
field description of h e f / ~ Within the mean-field model, 
h e ff is proportional to mean Mn spin-polarization (S) 
whose temperature dependence is given by the Brillouin 
function with a temperature dependent mean-field^ The 
curves in Fig. |2 can therefore be approximately assigned 
also to a (Ga,Mn)As DMS system where the effective field 
value changes through the temperature-dependent aver- 
age Mn-spin polarization rather than through the Mn- 
doping parameter x. Large values of h e ff correspond 
to low temperatures in this picture and, as seen from 
Fig. [3 the temperature-dependence of a in this regime is 
quite complex and sensitive to details of the DMS sample 
structure. Generally Fig. [21 suggests an initial increase 
of a with increasing temperature in samples with a large 
density of Mn moments, a nearly constant a for inter- 
mediate doping, and a suppression of a in samples with 
low Mn content. At high temperatures (small h e ff), the 
kinetic-exchange-driven Gilbert damping rate will grad- 
ually decrease towards zero with increasing temperature. 



FIG. 2: Gilbert damping coefficient a as a function of car- 
rier density for x = 2% - 8% and for quasiparticle life-time 
broadening of 150 meV (a) and 50 meV (b). 



III. EXPERIMENTAL FMR LINEWIDTH 

We now discuss our experimental FMR linewidth data 
recorded for the 120 nm thin Gai^^Mn^As layer with 
x = 8 % grown on GaAs (001) substrate. The FMR 
measurements were carried out at 9.46 GHz, with the ex- 
ternal dc magnetic field applied at different angles 9 with 
respect to the growth axis (9 = corresponds to the [001] 
crystal direction and 9 = 90° to the [110] direction). The 
Mn concentration in the sample was estimated from x- 
ray diffraction measurement of the lattice constant. The 
critical temperature in the as-grown (T c = 65 K) and 
annealed (T c = 110 K) samples were determined from 
SQUID magnetization measurements. A more detailed 
description of the sample properties and of our experi- 
mental set-up can be found elsewhere^ To analyze the 
measured peak-to-peak FMR linewidth AH pp , plotted in 
Fig. we recall the following general relation between 
AH pp and affl. 

AH pp (u;) = AH pp (0) + 4r— a ■ (11) 
v 3 gfJ-B 

Here AH pp (0) describes broadening due to sample 
inhomogeneity which is assumed to be frequency 
independen1i22i24 but can depend on the dc field orienta- 
tion. The second term in Eq. (|llf> arises from the Gilbert 
damping term in the LLG equation which gives a 
contribution to the FMR linewidth which is linearly pro- 
portional to uj and independent of the static magnetic 
field direction, if the magnetization is aligned with the 
fieldSi (and the dependence of a on the magnetization 
orientation can be neglected, as mentioned in the previ- 
ous sectionSi). This condition is satisfied in our sample 
since the FMR resonance field is larger than the magnetic 
field at which saturation magnetizations for different field 
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FIG. 3: Experimental peak-to-peak FMR linewidth in 
as-grown (filled symbols) and annealed (open symbols) 
Gao.92Mno.08As samples measured as a function of temper- 
ature for [001] and [110] dc magnetic field orientations (main 
plot) and as a function of the field angle at 4 Kelvin (inset). 



orientations coincided 

The strong dependence of the FMR linewidth on the 
field orientation in the as-grown sample suggests that 
magnetic inhomogeneities in the ferromagnetic layer con- 
tribute significantly to the FMR broadening. As seen in 
both the main panel and the inset of Fig. the angle de- 
pendence of AHpp becomes weaker in the annealed sam- 
ple and also the overall magnitude of AH pp is conspic- 
uously reduced. This observation is consistent with the 
improved quality of the sample (as indicated e.g. by the 
enhanced T c ) and implies that the leading contribution 
to the FMR linewidth might in this case come from the 
homogeneous (Gilbert damping) broadening. The right 
y-axis in the main plot of Fig. represents the exper- 
imental Gilbert coefficient obtained from the measured 
AH pp and from Eq. (jTT|) . assuming AH pp (0) — 0. Ex- 
perimental low-temperature values of a in the annealed 
sample are around 0.03. As seen from Fig. [21 these values 
of the Gilbert coefficient can be fully explained by the p-d 
kinetic-exchange mechanism of the damping of magneti- 
zation precession. However, because the density of Mn 
ions and their distribution in the lattice as well as the 
density of holes are not precisely known in our sample, a 
fully quantitative comparison between theory and experi- 
ment is not possible. The results suggest that experimen- 
tal studies of the frequency-dependent FMR linewidth in 
high quality DMS samples would be very valuable for un- 
derstanding the complex behavior of the Gilbert damping 
coefficient predicted in the theoretical part of this paper, 
and in separating those effects that are arising from the 
inhomogeneity within the epitaxially grown thin films. 



IV. CURRENT INDUCED MAGNETIZATION 
SWITCHING 

Using theoretical values for the Gilbert coefficient and 
anisotropy energy^ we now estimate the critical cur- 
rent for the spin-transfer induced magnetization switch- 
ing in a ferromagnetic semiconductor multilayer struc- 
ture. In general, spin-polarized perpcndicular-to-plane 
currents in magnetic multilayer systems can transfer spin 
between magnetic layers and exert current-dependent 
torques. For a trilayer structure, arguments based on the 
conservation of the angular momentum suggest the fol- 
lowing form of the torques on the two magnetic layers^ 
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where Mj^j = Mj^j/Af and I s is the spin-polarized 
electric current. The sign of the torque depends on 
the sign of the current, so that magnetization vectors 
in the two magnetic layers can be aligned parallel or 
anti-parallel by current flowing in one or the opposite 
direction. In a spin-valve structure with one "hard" and 
one "soft" magnetic layer, switching occurs when the 
torque in the soft layer overcomes the damping and the 
anisotropy terms. 

There have been a series of theoretical papers 1 - aimed 
at the quantitative description of spin currents and their 
effects on magnetization switching in metallic spin-valve 
structures. The theories are based on a two-channel 
model (spin-up and spin-down) and account for spin 
accumulation effects in the magnetic multilayers and 
spin transfer effects due to reflection at the ferromag- 
net /normal layer interface and due to the averaging 
mechanism associated with rapid precession of electron 
spins after entering the ferromagnet. The two-channel 
model is not applicable for semiconductor valence bands 
with strong spin-orbit coupling, complicating the quan- 
titative description of spin currents in DMSs. Strong 
spin-orbit coupling leads to a reduced spin-coherence 
time. However, the exchange coupling between the Mn- 
moments and hole carriers will make this time larger in 
ferromagnetic than in non-magnetic p-type semiconduc- 
tors. Experimentally, magnetic information can be trans- 
ported by charge carriers in DMS multilayers, despite 
strong spin-orbit coupling, as demonstrated, e.g., by the 
observation of the giant magnetoresistance effects^ 

For an order-of-magnitude estimate of the switching 
current in a (Ga,Mn)As-based magnetic trilayer struc- 
ture, we approximate the spin-current as I s ~ I( s ), 
where / is the electric current and (s) is the mean-field 
spin polarization of the itinerant holes in the (Ga,Mn)As 
layersiiS Adding the torque term (|12fl to Eq. Q for the 
soft magnetic layer we obtain an effective damping rate 
D = (Ua - I) /{I + a 2 ), where I = I{s)/(eN Mn SV). 
An instability occurs at D = and the corresponding 
critical current density for the magnetization switching 
is then given by j c = eUad/h(s) . Assuming a thickness 
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d oj 10 nm for the soft ferromagnetic layer and typical pa- 
rameters of a Gao.95Mng.05As DMS, U ~ 1 k J m~ 3 , a w 
0.02, and (s) w 0.3, the critical current j c ~ 10 5 A cm -2 . 
This estimate is two orders of magnitude smaller than 
critical switching currents characteristic of metallic spin- 
valves, 16 ' 17 primarily due to smaller saturation moments 
and anisotropy energies in the DMSs. Since the resistiv- 
ities are only 2 to 3 orders of magnitude larger in DMSs 
than in metals, observation of this effect should be exper- 
imentally feasible in a ferromagnet/normal/fcrromagnet 
semiconductor spin- valve structure^ 6 - 

We note that achieving low critical currents is par- 
ticularly important for magnetic tunnel junctions that 
are used in non- volatile memory devices. To avoid self- 
field effects, that lead to a spin vortex state rather than 
to a complete reversal and therefore to a smaller giant 
magneto- resistance effect^ the in-plane diameter r of 
metallic spin-valve devices with critical current densities 
j ~ 10 7 A cm" 2 must be of order ~ 100 nm>i& Mag- 
netic tunnel junctions have resistances too high for ap- 
plications when patterned to such small sizes. Since the 
Oersted field scales as ~ rj and the critical currents for 
spin-torque-induced switching we predict are two orders 
of magnitude smaller in DMSs than in metals, 1-10 /im 
size semiconductor tunnel junctions might still show suf- 
ficiently weak self-field effects and, therefore, a complete 
current-induced switching. 

V. SUMMARY 

In this article we have studied magnetization 
precession damping in ferromagnetic semiconductor 
(Ga,Mn)As alloys. We have attempted to employ the- 
oretical analysis combined with existing experimental in- 
formation to obtain the Gilbert damping coefficient and 
to predict the scale of critical currents for spin-transfer 
magnetization-reversal in these systems. In spin-transfer 
induced reversal, damping of magnetization precession 
competes with current-induced spin torques and deter- 
mines the scale of the current required to achieve rever- 
sal. 



Our theoretical analysis examines the mechanism that 
we expect to dominate at low temperatures in high- 
quality samples, due to the coupling of the d-level local 
moments to valence band holes. Wc derive an explicit 
expression for the Gilbert coefficient conventionally used 
to characterize damping in experimental studies, by com- 
paring microscopic linear response theory with the linear 
response limit of the phenomenological Landau-Lifshitz- 
Gilbert equations, and study how the values predicted 
by this model for the Gilbert damping coefficient depend 
on the hole density and on the size of the mean-field 
exchange interaction experienced by valence-band holes 
in the ferromagnetic state. We find that the magnitude 
predicted for the Gilbert coefficient, ~ 0.03, is consistent 
with experiment, but that the observed dependence on 
the external field and magnetization orientation is larger 
than can be accounted for by this mechanism. The ex- 
perimental FMR linewidth appears to have an inhomo- 
geneous broadening contribution that is not included in 
our theoretical modeling developed for homogeneous bulk 
systems. The uncertainty that presently exists in the 
relative importance of these two broadening mechanisms 
could be reduced by frequency-dependent FMR studies. 

In our view, the portion of the FMR linewidth broad- 
ening that is due to inhomogeneity should not, to a first 
approximation, be included in assessing the competition 
between spin-torques and spin-precession damping. We 
have therefore used our theoretical value for damping in a 
homogeneous system to estimate the critical currents re- 
quired for achieving magnetization reversal and obtained 
j c ~ 10 5 Acm~ 2 . 
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